Evidence of deep water penetration in silica during stress corrosion fracture 
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We measure the thickness of the heavy water layer trapped under the stress corrosion fracture surface of silica 
using neutron reflectivity experiments. We show that the penetration depth is 65-85 angstroms, suggesting the 
presence of a damaged zone of ~100 angstroms extending ahead of the crack tip during its propagation. This 
estimate of the size of the damaged zone is compatible with other recent results. 

PACS numbers: 61.05.fj, 81.40.Np, 68.35.bj 



If a small stone hits your windshield and you see cracks 
growing slowly from the impact, you are probably observing a 
stress corrosion process. Under very moderate external tensile 
stresses, cracks may indeed grow in silicate glasses, thanks 
to a chemical reaction which involves the water molecules of 
the surrounding environment. This is a complex phenomenon, 
which started to be studied in the sixties lH] |2l and is not yet 
fully understood (see [T| for a recent review). As a matter of 
fact, two important questions remain to be solved: the exact 
mechanism by which water molecules manage to break Si-O 
bonds, and the role of the amorphous structure of glass in the 
fracture properties. 

In the classical picture, proposed first by Michalske and 
Bunker ID, water molecules break the Si-O bonds located 
exactly at the crack tip thanks to a hydrolysis reaction. For 
small enough external loads, the crack velocity is controlled 
by the rate of the chemical reaction, which depends both on 
the degree of ambient humidity and on the applied stress. This 
regime is traditionally refered to as "Stage I" |[l]|5l. At higher 
applied loads, when the crack velocity reaches a value that de- 
pends on humidity, the slowest phenomenon (which imposes 
its kinetics to the crack velocity) is the diffusion of water 
molecules to the tip along the fracture surfaces. Since sur- 
face diffusion is not sensitive to the external applied load, the 
crack velocity in this "Stage 11" does not depend on it either. 

This classical picture, however, does not take into account 
the disordered character of the glass structure. In a perfect 
crystal, where atomic bond orientations and energies are (5- 
distributed, bonds at the crack tip will break first because 
stress concentration is maximum there. But when consider- 
ing a more complex arrangement of chemical bonds, disjunc- 
tions are likely to occur at some distance away from the tip. 
This effect was seen first in Molecular Dynamics (MD) simu- 
lations of dynamic fracture, where a nanometric Process Zone 
(PZ) was observed to form ahead of the main crack front |6 
|8|. Somewhat later, such a "quasi-brittle" ll5]|9l behavior was 
claimed to be observed experimentally using in situ Atomic 
Force Microscopy (AFM) experiments lfT0l4T2ll . These obser- 
vations are however still very controversial llT3l[T4ll . Because 



AFM observations are restricted to the free surface of the 
specimens, several artifacts can alter the measurements [i5|. 
Moreover, results may be tampered by tip effects |16|. As a 
matter of fact, there are several significant differences between 
the free surface and the bulk of the specimen. Of particular 
importance is the exposure to water in the case of stress cor- 
rosion fracture: while the free surface is in direct contact with 
the ambient humidity, water molecules have to diffuse within 
the material for Si-O bonds to break at a distance from the 
crack tip. Although experiments have been performed at high 
temperature only ifTTlfTSl . a rough extrapolation of Tomozawa 
et al's results 1,1^ suggests a water diffusion coefficient in 
silica of the order of ^lO^^^cm^.s^^ at room temperature. 
This means that the penetration length of water molecules 
into unstrained glass should be approximately 3pm (respec- 
tively 0.3 A) during the time it takes for a crack moving at 
lO^^m.s^^ (respectively at 10~^m.s^^) to cover 100/im. 

However, because of the huge stresses concentrated at the 
crack tip, diffusion is enhanced by orders of magnitude in the 
vicinity of the tip during fracture 112011211 . as observed in sev- 
eral other materials [22.. ,23 1. Therefore, water is expected to 
penetrate into the glass and, because of the heteogeneity of 
the material mentioned above, start breaking bonds and cre- 
ate microcracks ahead of the crack tip. This in turn increases 
further the diffusion of water, thereby creating more corro- 
sion and potentially leading to a substantial damaged zone. If 
this scenario is correct, a rather thick layer of water should 
remain trapped underneath the nominal fracture surface after 
the crack has propagated and stresses have relaxed. Since the 
diffusion constant is so small in unloaded silica glass (more 
than 100 days for travelling 1 nm), one should observe post- 
mortem a "fossil" water profile, essentially frozen-in at the 
time of its creation, with a thickness of the order of the size of 
the damaged zone. The aim of this work is to provide quan- 
titative evidence for the above scenario using neutron reflec- 
tivity 1 24, 25 1 to measure the thickness of the water layer left 
behind the crack. We find that the penetration depth of water 
is of the order of a hundred angstroms, suggesting the pres- 
ence of a large damaged zone. 
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FIG. 1. (a) Picture of the experimental setup. The neutron beam 
is schematized in red, incoming from the last slit of the collimator, 
reflecting on the sample and going within the slit in front of the de- 
tector, (b) Picture of the broken sample, showing the two black sheets 
of B4C used to select one area of interest (see text), (c) Sketch of the 
broken sample, with the two parts described in the text. 



Fracture experiments. Fracture experiments were con- 
ducted in a highly controlled manner via Double Cleavage 
Drilled Compression (DCDC) samples. In this geometry, the 
stress at the crack tip naturally decreases, enabling us to con- 
duct all our experiments in the stress corrosion regime. DCDC 
samples used herein are cuboids of size 5 x 5 x 25mm'^ with 
a 1mm diameter hole drilled in the center. They are made 
of Corning 7980 pure silica. The fracture experiments were 
conducted in a glove box which had been saturated with 
heavy water. Two symmetrical precracks are first initiated 
from the hole of the sample as described in ifTTI . Subse- 
quently the load is adjusted in order to reach a desired ve- 
locity |11|. Zone 1 (Fig. [T];) corresponds to a stress inten- 
sity factor 0.61MPa.m^/-^. Zone 2 (Fig. [l];) corresponds to 
a stress intensity factor 0.77MPa.m^/^ (average velocity of 
4.10-^m.s-i). 

Heavy water has been chosen because its coherent length 
density — 6.39 lO^^A^^ is higher than the one of silica 
{hs = 3.41 10~^A~^), and of opposite sign to that of light 
water (—0.53 IG^^A^^). If some water is trapped in the vicin- 
ity of the surface of the sample, the reflectivity of the sample 
should thus increase in the presence of heavy water whereas it 
would only weakly decrease with light water. 

Neutron reflectivity experiments. Specular Neutron Re- 
flectivity (SNR) measurements have been carried out on the 
horizontal time-of-flight EROS reflectometer (Saclay, France) 
with a fixed angle of 1.195, with a neutron white beam cov- 
ering wavelengths A from 4A to 25 A, covering an accessible 
q-range (diffusion vector q = 27rsin0/A) from 0.005 A^^ to 
O.O32A-1. 

Zones 1 and 2 (see Fig.[T];)) described above were studied. 
In order to select one of these areas of interest, we used the 
following trick. The sample was almost completely hidden on 
the neutrons path by two black sheets of B4C, a strong neutron 
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FIG. 2. Experimental neutron reflectivities plotted as a function of 
the diffusion vector q for: the control sample (*), the specimen bro- 
ken at a stress intensity factor Ki — 0.61MPa.m^^^ (Zone 1, * 
) and for the specimen broken at a stress intensity factor Ki = 
0.77MPa.m^^'^ (Zone 2, + ). We have superimposed to the ex- 
perimental curve corresponding to the control sample the theoret- 
ical Fresnel reflectivity curve (red dotted line) as well as the re- 
sult of the second order Born approximation (blue dashed dotted 
line). The experimental results corresponding to Zone 1 and Zone 
2 are fitted using Eqs. l[3](, l|5j and l[6](. The best fits correspond to 
00 = 0.348 ± 0.003, T ^ 43A and « 35A in Zone I, and 
= 0.567 ± 0.003, £" ^ 46A and A" ^ 23A in Zone 2. Inset: 
Contour lines of the fit root-mean square error in the plane ((., A) in 
Zone 2, showing that while the combination £ef / = A + £ is rather 
well pinned down by the fit, A — ^ is a "soft" direction. The rela- 
tive experimental rms error per point is 0.075, whereas the minimum 
relative error achieved by the fit is 0.06. 



absorber, to let the neutrons illuminate only the desired region 
(Fig.[T]a and b). In order to test that the selected region was flat 
enough to allow a correct measurement, we have checked that 
the half full width of the alignment rocking curve was lower 
than 0.25°. When this was not the case, the illuminated region 
was reduced until this condition was met. The resulting illu- 
minated surfaces were very small, of the order of ^ 25mm^. 
Because of this smallness, we used very long counting times 
to get a reasonable noise-to-signal ratio, up to 48 hours per il- 
luminated region. In particular, we measured the background 
independently from the sample with great precision, enabling 
its subtraction with a good accuracy. 

Results. The experimental curves presented in Fig.[2]clearly 
show a huge change in the reflectivity of broken samples when 
compared to the reflectivity of an unbroken control specimen. 
Since the roughness of the fracture surfaces is larger than the 
roughness of the control specimen, one should a priori ex- 
pect a decrease of the reflectivity of the broken samples. The 
difference seen in Fig.[2]is therefore underestimated. 

Fig. |2] shows also that the reflectivity of the control spec- 
imen corresponds perfectly to the Fresnel reflectivity Rp = 
{q — {q^ — q^))"^ / {q + \/ (g^ — <?c))^ of a semi-infinite sil- 
ica diopter for which the coherent length density is equal to 
6s 1251 l26ll . for which q^ = \/4ttFs. In order to fit the neu- 
tron reflectivities measured in Zones 1 and 2 (Fig. |2]l, we 
have used a second order Bom approximation, assuming that 
the heavy water concentration is not a constant, but that it 
decreases with the distance z from the free surface of the 
tested sample as cj){z). This concentration profile translates 
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into a coherent length density profile bw{z) = hw4){z). Since 
most of the workable signal is obtained in a region of large 
diffusion vectors q, far enough from total reflexion, reflec- 
tivities are quite small, and hence the reference situation is 
the free case, when all the incident neutrons are transmit- 
ted. The neutron wave function ?/; obeys the following eigen- 
value equation: (Pip/dz'^ — —q^ip + V{z)ip, with V{z) — 
47r(6s + b^(p{z)) a small perturbation: V(z) ^ q^. We write: 
4> — ipo + i'l + '4'2 + with il^i and the first and second 
order corrections (in V). Adapting the calculation of ll27l . 
one finds: %pn+i{z) = dz'V{z')ipn{z')Go{z, z'), where 
the Green function Go{z, z' > 0) relevant for our boundary 
conditions is: 

Go{z,z') = -^cxp {-iq{z' - z)) ifz<0 
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exp {—iq\z' — z\) ifz>0 



(1) 



Using the above expressions, we get the expression of the re- 
flectance r to order V'^: 
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The reflectivity Ris R — | ?- ^ . In order to check the validity 
of our second order Born approximation, we first verify that 
our result in Eq. Q tends to the Fresnel reflectivity for high 
values of q when is a constant equal to q^ = A-rrbs- This 
limit leads to the following reflectivity: 
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which coincides with the corresponding large q expansion of 
the Fresnel reflectivity (see Fig.|2]l. 

We then tried to fit the reflectivities in Zones 1 and 2 us- 
ing the simplest function involving a single length scale, i.e. 
(j){z) — (po exp (— z/A). Although this can be made to fit the 
Zone 1 results, the reflectivity increase in Zone 2 is too large 
to be accounted for using this simple function. Hence, guided 
by the idea that there might be a saturated layer of depth £ 
close to the surface that becomes more diffuse deeper in the 
sample, we posit that: 



4>wiz) = 00 if z < ^ 

0^(z) = 0oexp(-(z-^)/A) ifz>^ 



(4) 



This choice leads to a reflectance r that can be written as: 
r = ro + ri + r2, with tq as the Fresnel reflectance (Eq. Q) 
and: 
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Fig. |2] shows the best fits of the experimental measure- 
ments performed on the two fracture surfaces using R ~ 
\ro + ri + (Eqs. (jsjl, (jsjl and The best fit is achieved 

' =rf 43A and A^ « 35A in 



with 



= 0.348 ± 0.003, 



(slow) Zone 1, while in (fast) Zone 2, 0" = 0.567 ± 0.003, 



« 46A and A^^ « 23A. Note that (f>Q is veiy accurately 
determined by the fit, although the error bar we quote only 
accounts for statistical uncertainty, and not systematic effects 
coming from the choice of the fitting function and of the inter- 
val over which the data is fitted. On the other hand, the quality 
of fit has a "soft direction" in the plane (£, A), as represented 
in the inset of Fig. [2] As expected, the total effective width 
of the layer, £ef / = £ + Ais better determined than £ and A 
separately. The statistical error bar on £ef / is smaller than 1 A, 
but again systematic errors are much larger. As shown in the 
inset of Fig. |2] £ef / can be varied by ^ ±5 A and still lead to 
an acceptable fit. 

Discussion. Our results therefore clearly show that heavy 
water is present over ~ 65 to 85A under the stress corrosion 
fracture surface of pure silica in the conditions of our exper- 
iments. This penetration depth is much larger than what is 
expected using a room temperature extrapolation of the diffu- 
sion coefficient of light water in silica |,19J (to our knowledge, 
there are no such results concerning heavy water, for which 
diffusion should be even smaller). A likely explanation is as 
follows: diffusion enhancement in the vicinity of the crack tip, 
where huge tensile stresses are present ( ll20l - ||231 ). allows 
the water to penetrate inside the bulk and create a damaged 
zone which helps water progressing further still. Our neu- 
tron scattering experiment gives information about the water 
penetration depth in the direction perpendicular to the fracture 
surface. To investigate how far water penetrates parallel to the 
crack propagation direction, one would need a detailed self- 
consistent model for the coupled growth of the damaged zone 
and the diffusion of water. We grossly simplify the problem by 
computing the dilation field D{x, y) — e^xix, y) + £yy{x, y) 
in the vicinity of the crack tip within a purely elastic model, 
and postulating that the water penetrates in a region defined 
by D{x, y) > Dc, and infer the anisotropy of the water pene- 
tration from that of the iso-dilation lines. We analyze DCDC 
samples identical to those used in the experiments and com- 
pute the dilation field in the mid-plane of the specimen us- 
ing finite element simulations, with element sizes decreasing 
exponentially as approaching the crack tip, so that the strain 
field is resolved at the nanometer scale within that region. As 
shown in Fig. 3 for an external loading K = 0.61 MPa.m^/^, 
the domain at the crack tip with a high level of dilation ex- 
tends deeper in the direction of propagation than in the per- 
pedicular direction, by a factor y/2. Therefore, we estimate 
that water penetrates roughly 9 — 12nm ahead of the crack tip, 
in the crack propagation direction. If one identifies the water- 
rich region with a damaged zone, our estimate is in agreement 
with the fact that the strain field observed in the vicinity of a 
stress corrosion crack tip is elastic only on scales larger than 
- lOnm |28l. 

We observe an increase of the neutron reflectivity, not only 
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FIG. 3. Finite element simulations of the dilation field D{x, y) = 
^xx {x,y) + eyy {x, y) in the vicinity of the crack tip in the mid-plane 
of the specimen for an external loading K = 0.61 MPa.m^'^'^. The 
crack propagates from right to left, and the color code indicates the 
amplitude of D(a;, y). 



with respect to the Fresnel case, but also with respect to the 
case where we assume a heavy water single-length scale, ex- 
ponentially decreasing profile. This has prompted us to postu- 
late the presence of an homogeneous layer of water, followed 
by an exponentially decaying profile. The width of the layer, 
i ss 4nm, is significantly larger than the height fluctuations on 
silica fracture surfaces which do not exceed Inm [29, 30 1. It 
is tempting to interpret this zone as a strongly damaged zone, 
with a density of microcracks that is larger when the stress in- 
tensity factor - and the crack velocity - is higher. This would 
naturally explain why > 0q, although we have no clear 
explanation as why w . This might be due to the cancel- 
lation of two opposite effects: a stronger stress enhances the 
diffusion of water in the bulk, but at the same time the crack 
speed is larger, leaving less time for the corrosion mechanism 
to operate. The latter mechanism in fact explains why the 
exponential region, which is probably more sensitive to dif- 
fusion, is wider in Zone 1 than in Zone 2 (A^ > A^^). In 
order to be more quantitative, one requires at this stage a de- 
tailed model for the formation of the damaged zone, dynam- 
ically coupled to the water profile. A more systematic study 
of the effect of the external applied stress on the heavy wa- 
ter content would be needed to provide a sound basis for such 
a quantitative model. Because the damaged zone size is ex- 
pected to be even larger for lower crack velocities |29 3^, 
further experiments will be performed for crack speeds as low 
as 10~^^ — 10~^°m.s^^, with larger specimens such as to in- 
crease the intensity of the reflected neutron beam. 
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